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SUMMARY 

A higher order theory of the motion of close satellites in the vicin- 
ity of the critical inclination is developed for the case of small ec- 
centricity as well as moderate eccentricity. Terms up to the third 
order of magnitude are included for a potential which is plane sym- 
metrical. It is shown that solutions up to this order can be obtained by 
using only Jacobian elliptic functions for the case of small eccentricity. 
However, for moderate eccentricity, elliptic integrals of the second 
kind or the Jacobian zeta- function are required. To obtain further ap- 
proximations, elliptic integrals of the third kind are necessary for the 
case of small eccentricity; but, for the case of moderate eccentricity 
only the Jacobian elliptic functions and elliptic integrals of the second 
kind are required. In any case, the solutions are, in general, periodic 
functions— as is well known from the form of the equations of motion. 
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A CONTRIBUTION TO THE THEORY OF CRITICAL 
INCLINATION OF CLOSE EARTH SATELLITES* 

by 

Shinko Aoki 1 

Goddard Space Flight Center 

INTRODUCTION 

Many authors have shown their interest in the motion of close earth satellites near the critical 
inclination (for example, References 1-5). However, almost all of these authors have considered only 
the first terms, which are of the order of the square of j 2 in the equations of motion. However, it was 
first shown by Izsak (Reference 5) that continuations of the ordinary treatment of the libration to the 
higher approximations would break down; and, in the case of small eccentricity, some libration of a 
peculiar kind including the next order of magnitude would occur. The first term, which has e 0 2 as a 
factor, is so small that the next order term— which would be considered as a higher order of magni- 
tude in the normal case— becomes the same order of magnitude as the preceding term. However, he 
presented this only for a consideration of the form of the Hamiltonian; he did not try to solve the 
equations of motion with time as the independent variable. 

Here some expressions for these solutions will be presented, not only for the case of small ec- 
centricity, but also for the normal case including terms up to the third order of J 2 in the original 
equations of motion, where j 2 is the coefficient of second zonal harmonics of the earth's potential and 
j 4 is assumed to be of the order of the square of j 2 and, further, J 6 to be of the third order of magni- 
tude. The odd harmonics will be totally neglected. A higher order theory including these terms is 
under consideration, but must be treated in the future. 

The method adopted is to introduce intermediary solutions, of which the argument \u is only 
slightly different from time (the independent variable) except for a constant factor, and such that by 
which the amplitude of the libration or the change of states— in other words, a transition from revolu- 
tion to libration— is fully determined. In the case of small eccentricity, the intermediary solution can 
be obtained by the following integration: 

. = \ du , 

y(l - a sin 2 (l - sin 2 oij) 


•Published in substantially the same form in Astronom. J. 68(7): 35 5-364, August 1963. 
fNAS-NASA Research Associate. 
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where n/2-co 1 is the mean value of the argument of perigee and K is constant; this integration can be 
performed even if only the Jacobian elliptic functions are used, as will be seen in a later section.* 


Further approximations, beyond the aforementioned order of magnitude, would require more com- 
plicated calculations. Some comments will be made on the subject although explicit formulas are not 
included. Izsak, after removing the so-called periodic perturbations, wrote the equations of motion: 


dp 


dr 

dco 


d co 


dr 

dp 


where <u is the mean value of the argument of perigee; p is some function of G, the conjugate variable 
to co; and 


<J> - A(p) + j B(p) cos 2 co , 

A (p) = p 2 + 2 a n p n , 

B Ip) = Z£ n P n . 

0 


/3 0 = e 2 [l+0(j>] , 


/Sj = - 2 [l + 0( j )] , 


l/J 4 

i = -lljT + J 2 j Po" > 0 - 


3 J 2 n o 
CPo 2 


Since r has a factor of J 2 , the terms of the order of j in o> correspond to those of the order of J 2 2 or 
j 4 in the original equations of motion. 

In this paper we have started with a Hamiltonian similar to the one used by Izsak but have added 
some additional terms in order to more easily consider the higher order terms; although this explicit 
evaluation, from the coefficients of zonal harmonics of the earth's potential and from the integration 


*A list of symbols is given in Appendix A. 
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constants, is not given. We are concerned with the solutions of the equations of motion which contain 
the Hamiltonian up to the order of j 2 : 

<E> = A( / o)+jB( ( o)cos 2 <i>+j 2 C(p) cos 4 co , 

although some comments will be made on obtaining higher order approximations. 


EQUATIONS OF MOTION 

We shall now consider the following equations: 


where 


d£ ... 


dr 

dco ’ 

d co 


dr 

~ dp 


$ = 


A ip) 

B (p) 
C(p) 


A ip) + j B(,o) cos 2 a> + j 2 C(p) cos 4 &j 

CO 

E a - pn ’ 


n=2 


= L p P n 
0 n 


= Ir n /> n 
0 " 


(1) 


and a n , P n , y n are constants which satisfy a 2 > 0, / 3 0 > 0. 

In order to remove the first degree term of p in the Hamiltonian function, we change the variable 
p to x = P - p (cos 2 «) such that 


d x)x=0 


d<$>\ 

dp) 


= 0 


P=p 


( 2 ) 


If we take the terms up to j 2 in Equation 2, we have 


j /5 1 cos 2 oj + j 2 7 j 


+ 2 


(a, + j fi cos 2 


" P 


3 a 3 P 2 = 


3 



from which 


Thus, the equations, 


1 fii 


. - . 1 Vl l ?! $2 3 a 3 P\ i . 

1 I — ~ 4- _ «— y- I r'rvp ‘ , 


2 J — c °s"o) + 2 — + 2 , “ 8 ’ icos’o. 


(3) 


dx 

dr 


d^i 

dr 


dco ' 


dv 

dy 


(4) 


are obtained with 


W( x , oj) = 4>(p+X. ") 


3 0 + a 2 X 2 + a 3 X 3 + a 4 X 4 


where 


a o - i /3 q cos2 “ + j 2 (" T ^l 2 a 2 _1 + ?o) cos4 w + 


a„ - a. 


+ j l/3 2 - -J a 3 Z 3 ! a 2 ' j cos “ 


(5) 


a 3 + 


a 4 = <x 4 + 


We have considered terms up to the order of j 2 in >P, assuming x is the order of j 1/2 . We know 

that 


v = jr , 


(6) 


where r is constant, is an integral of Equation 4. Neglecting higher order terms, we can easily obtain 
an expression for x from Equation 6: 


X 


Xj 


l 

- 2 a 3 


a 


-l 

2 


Xi 


+ 


4-0 o-i 



Xi 


(V 
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where 


(i r_a o) a 2 _1 


Rewriting x, 2 


where 


and putting co t 


2a 2 x 


where 


Xi 


a, j 


r-/3 0 cos 2 CO + j ^r /3 2 a 

+ j (r-,fl 0 cos 2 co) (- / 3 2 <x 2 _ 1 + ja 3 / 3 , a 2 ' 2 j cos 2 a> 


-i 1 -y 0 ) cos4 “ 

3 


in the form 


Xj 2 = fjL ^1 + Pj COS 2 CO + Qj COS 4 , 


~ r a 2 _1 3 - 


r -l 


/ 3 0 + j r (“^2 a' + 2 ' a 3 ^i a 2 


Qj j r 1 (4 fil a 2 1 To + ^0 a 2 2 a 3 ^0 a 2 


77/2 - co, we have 

dco. 


dco. 


+ 3a 3 y 2 + 4a 4 x 3 2a,ii W 


^ X 0 


1 + a 3 a 2 1 A 1 1/2 Xo + (2 a 4 a 2 1 “ 8 a 3 a 2 ) 


^Xo 


( 8 ) 


dr , 


— - 1 - 1/2 

Xo ~ a 2 a 2 A 1 Xi 

= ^1 + Psin 2 a)j +Qsin 4 cOj 

P “ fio ^ 1 + j(®2 a 2 2 a 3 a 2 ^ 

Q = j r 1 (4" Z 3 ! 2 a 2 1 ""Xo ” ^0 <®2 a 2 1 + T a 3 ^0 a 2 ) • 


> 


J 


0) 
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INTERMEDIARY SOLUTIONS 

From the denominator on the right side of Equation 8 we can see that it will be zero if, and only 
Xo ~ Provided ^ = 0( j r ) is small. This means that the essential character of the trajectory is 
determined entirely by x 0 . In the case when a> 1 is real— which includes, in general, j 2 terms of the 
equation of motion. Thus, Xo discriminates among three kinds of motion: libration, revolution, and 
the critical case. 

Therefore, we first must consider 


dan 


|/l + P sin 2 CiJj + Qsin 4 Wj 


Putting sin 2 ojj = x , we have 


do)j 

y (l “ a sin 2 CL> X j (1 - fi sin aij 
du . 


( 10 ) 


du - 


dx 

2 |/( 1 - x) x (l +Px + Qx 2 ) 



and, putting y = x _1 , we have 


du - 


~ dy 

}/4!y- 1) (y 2 + Py + q) 


~ dy 

V4(y- 1) ty-a) !y-/3) ’ 


or 


du 


~ dz 

]/ 4 ( z " e f) ( z_e g ) ( 2 ~ e h ) 


where 


( 11 ) 


and 


P- l 

z = y + — 3 - , 


Thus, we have 


e , e 
g 


0 = 


+ T- 


a + 


P- 1 


p + 


p - 1 


z 
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*°( u-u o) , 



where #>u is the Weierstrassian elliptic function of the second order. We assume that the integration 
constant is such that 


U 


0 


0 , 


which corresponds to the initial condition of a>j = 0 when u = 0, as will be seen later. Then we have 

/ p - i 1 

sin 2 co 1 = Uu 3—j = (#u - e f + l)" 1 

= (l + fj 2 u) _1 , 
or 


cot Ci) J - ± f j u , 


and 


yfl + P sin 2 ojj + Q sin 4 co 1 - gj u • hj u (l + f j 2 u) " 1 . 

In the above expressions, f ju is an associated primitive elliptic function of the first order (Reference 
6), whose poles are u = 0 and 2&> g and whose zeros are u = a { and co { + 2u g within a parallelogram 
(2 co f , 4w g ) , where we have 

/%CO 

dz 

Je -j/4 (z - e f ) (z-e g ) (z - e h ) 

P 

with p = f , g, h and co { + + cj h =0. The other gj u and hj u are given by permutations of (f , g , h). 

However, this expression is so general that, at first glance, we cannot see hpw the function be- 
haves for real values of u, which is the only case of interest. Therefore, in order to clarify its 
behavior, we shall divide the problem into several cases. 

(i) a , j3 real: Libration and critical cases 

(a) 1 < /3 < a < y , libration within - a ~ 1/2 < sinw, < a - 1/2 

Since 


pu 


1 


+ 


P-1 


a - 1 


3 


+ 


sn 2 (\.u, k) 


( 12 ) 



where A. 2 = a- 1 and k 2 = ip- l)/(a- 1), we have 


(\u, k) 


(\u, k) 


and 


( 13 ) 


Xo 


y 1 + Psin 2 *)j + Qsin 4 ^j 


jg - 1) cnXudnAu 
a ~ cn 2 (A. u, k) 


The case 1 <y </3<a does not correspond to any real case in the theory of the critical inclina- 
tion of close artificial satellites (Reference 5); therefore, we omit this case here, although similar 
expressions are obtainable if we substitute sinu = cos co t instead of sin co 1 in Equation 1. 

(b) p< l < a< y , libration within - a " 1/2 < sin co t < a ~ 1/2 

Since 


fpu 


F - 1 a ~ /3 

sn 2 (Xu, k) 


where X 2 = a- p andk 2 = (l -p)/[ a -p), we have 


sin 2 a 


sn 2 (Xu, k) 
a - /S cn 2 (Xu, k) 


and 


{a- /3) cn ( Xu , k ) 
a -/Sen 2 (Xu, k) 

(C) p < a < 1 < y, revolution 


Since 


£>U 


P + 


P- 1 


1-/3 


(Xu, k) 


where X 2 = l-^andk 2 = (a -p)/[ 1 - p), we have 


sin 2 co j 


sn 2 (Xu, k) 

1 ~Pcn 2 (Xu, k) 
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and 


( 1 - f3 ) dn ( Xu , k ) 
1 - /3 cn 2 ( Xu , k ) 


Xo “ 

(a') l</3 = a < y, critical 
Since sn\u - tanhXu when k - 1, we have 


sin 2 ojj 


sinh 2 Xu 
a cosh 2 Xu - 1 


where X 2 = a - l and 


a ~ 1 

a cosh 2 Xu - 1 

(b ') /3 < l - a < y, critical 

From either case b or c, when k - l, we have 


and 


sin 2 oij 


sin h 2 Xu 
cosh 2 Xu - (3 


1 1 - /3) cosh Xu 
cosh 2 Xu - /? 


where X 2 = 1-/3. 


(ii) a, ji imaginary: Revolution 

Here the direct transformation from g >- functions to Jacobian elliptic functions does not give a 
real function. Therefore, another type of transformation is used, as is described by Whittaker 
and Watson (Reference 7, p. 513), since p and Q are real. 

If we apply the transformation 


where 


y 


c 


x-B 
x-A ’ 


A 


- 1 - V 1 +P + Q 
P + Q 
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B 


- 1 + y i + p + q 

P + Q 


„ _ - 2 -P-Q -2 yr+p + o 

c P + Q 

to Equation 11, we then have 


ck), 

d U = — — — — 

+ Psin 2 o ; 1 + Qsin 4 cj 1 

= dy 

2X -/(l-y 2 ) (k' 2 +k 2 y 2 ) 

where 

^ 2 = yi + p + q = y ( i - a) d-/3) , 

u2 = -P-2 + 2 yi+P+Q _ a + fi - 2 , 1 

4 yi + P+Q 4 y(l-a) (1-/3) + 2 

Accordingly, we obtain 


y - cn [2X. (u - u Q ) , k] . (14) 

It should be noted that the form of y in Equation 14 yields a real function, when \ 2 > o and 1 > k 2 > o . 
This is made possible by choosing an appropriate sign of the above root function (which is, of 
course, a double-valued function) in such a way that 


yi + P + Q > 0 ; 

this is possible because P and Q are real and satisfy the inequalities 

l V 

1 + P + Q = (l + -f) - 4 (P 2 - 4Q) > ° , 

(P + 2) 2 < 4( 1 + P + Q) , 


since 


A = p 2 - 4Q < 0 . 
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Table 1 


List of Solutions for daij - y g du , y g - ( 1 - ax ) ( 1 - /3x) , x - sin 2 - cos 2 co , 


Case 

Designation 

X 

Xo 

k 2 

X 2 

Character of the Motion 

ia 

1 £ P 5 “ 5 x' 1 

' 

sn 2 (Xu, k) 

(a - 1 ) cn Xu dn Xu 

P~ 1 

a - 1 

libration between 

- a~ 1/2 < sin&jj < a ~ 1/2 

a - cn 2 (Xu, k ) 

a - cn 2 (Xu, k ) 

a - 1 



sn 2 (Xu, k) 

( a - P ) cn Xu 

1-/3 

a - p 


ib 

P < 1 < a < x _1 

a - P cn 2 (Xu, k ) 

a - p cn 2 (Xu, k) 

a “ 

1 ibration between 

- a ~ 1/2 < sintdj < a _1/2 


ic 

P < a < 1 < x' 1 

sn 2 (Xu, k) 

( 1 - P ) dn Xu 

a - P 


revolution 

1 - /S on 2 (A.u, k ) 

: 

1-/3 cn 2 (Xu, k) 

1-/3 

1 P 

ia' 

1 < P = a < x' 1 

sinh 2 Xu 

a - 1 



critical, limiting case of /a. 


a cosh 2 Xu - 1 

a cosh 2 Xu - 1 


a— 1 

when k 2 - 1 

ib’ 

P < 1 = a < x _1 

sin h 2 Xu 

( 1 - P ) cos h Xu 


i - p 

critical, limiting case of both 


cosh 2 Xu - P 

cosh 2 Xu - P 


ib and ic , when k 2 -* 1 

ii 

a / P, imaginary 

sn 2 (Xu, k) 

X 2 (l - k 2 sn 2 Xu cd 2 Xu) 

1/ 2-a-p \ 

V( l-a)(l-/3) 

revolution 


1 < x" 1 

1 - cn 2 Xu + X 2 cd 2 Xu 

1 - cn 2 Xu + X 2 cd 2 Xu 

2 K 2 V ( l-a) ( \~P) j 


Hi 

a = /3 = 1 < x" 1 

u 2 

1 



critical 

1 + u 2 

1 + u 2 












THE CASE OF SMALL ECCENTRICITY 

If /3 0 (%e 0 2 ) and r are of the same order of magnitude as j ,then the argument given for the equa- 
tions of motion becomes somewhat simpler because x is of the order of j so that the terms a 3 x 3 and 
a 4 x 4 in Equations 5 are negligible provided we include terms up to the second order in these equations. 
Moreover, in the expression ^xo 2 it is sufficient to consider only terms up to the second order. In 
the case of small eccentricity near the critical inclination, y 0 itself (Equations 5) would be of the 
order of j 2 , which is also negligible (Reference 8). 

Thus, in this case, the intermediary solutions described in a previous section would give suf- 
ficient accuracy provided that we take Equations 1 to the second order. 

Now putting /3 0 = j/3' , r - jT' we have 

p - 
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Q - \ r' -1 /3j 2 a" 1 , 
/a = j 2 r'a 2 _1 : 

from which we obtain 

(a,/3) = (l ± 

From Equation 8 we have 


u - 2a, y^T r 


where 


H - j 2 r'a 2 -1 


If T' > /3 0 ' 2 a 2 /3^ 2 , a and /3 become imaginary so that this case corresponds to the Case ii . 

On the other hand,r < /3 0 ' 2 a 2 /3" 2 corresponds to Case i or Case Hi. The subdivision within Case i 
can be obtained easily 


if 


0, 2 P 0 ^ 

Po - 4^1 < r ' < T - then 


Case ia; 


0 < r ' < P 0 - 4T 2 - then 


Case ib; 


and if 


F ' > P 0 


P 2 . Po 


4a . 


and -y , then 


Case fc. 


Rewriting the above and combining it with the real/imaginary criterion results in the following scheme: 


Case A 


P? Po , P 2 Po 2 «2 

Po>^- 2 OV — < Po-^ 2 < -Jf 
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if 


Case ib; 


if 


if 


Case B 


/3 2 

0 < r ' < K - 4^7 ’ then 


/V 


IL 

4a. 


< r' 



1 then 


lo ,2 « 2 

P 2 


< 


P', then 


Case ic; 


Case ii. 


if 


if 


N 


1L 

2a 2 


or 




ll_ Pq 2 «2 _/v 

4 a 2 < K 2 


0 < r < /3 0 ' 


4^7 ’ then 


^0 


il 

4a, 


< r' <- 


Po 


/V 


• then 


Case ib; 


Case ia; 


Po 2 -2 , 

— TJ— < r , then Case ii. 

P-1 

The critical cases are easily obtained by comparing the above scheme with Table 1. In order to be 
able to distinguish between the different cases, a diagram is drawn with the parameters expressed by 
P' and /3 0 ' in units of /3,^/4a 2 (Figure 1). 

Now, the period of sinwj in terms of r is given by 


T 


4K 

2a" 2 j T ,1/2 \ ’ 


where the complete elliptic integral of the first kind, 


K 


/» 7T/2 


do)j 


J 0 


y 1 - k 2 sin 2 co l 



(15) 


Accordingly, if in Figure 1 we take a point which is moving around the critical point in in a 
counterclockwise direction, we will find the change of the period and of the character of the motion as 
given in Table 2. 


14 



3 


In any case the solutions of Equations 1 up 
to the order considered here, are given in case 
ia by using Equations 3, 7, 9, 12 and 13; 


sn 2 Xu 
a - cn 2 Xu 


p = p+ x - p+ M 1/2 Xq 


.1 • sn 2 Xu 

2 1 a 2 a - cn 2 Xu 


+ j 



( a - 1 ) cn Xu dn Xu 
a - cn 2 Xu 


where 



Figure 1— Criterion for several cases, expressed by pa- 
rameters T' and /3 0 ' in units of (l/4)^S| 2 a 2 " 1 j The respec- 
tive cases are described in Table 1 . 


a 



r '/W 



k 2 - ^ \ ’ Xu " 2 ia - 1 a 2 1/2 r' 1/2 jr . 

For the other cases, we can easily obtain the results from Table 1 in a similar manner. 

Table 2 

Behavior of the Period and the Character of the Motion. 


s) 


Case 

a 

a' 

a 

b 

b 

b' 

c 

c 

a 


/ 


X 

2rr 

y 



2n 

S' 

Period T 

/ 

CO 

X iGT- 

1) T' • 2a 2 

j / 

CO 

x vnr 

- a ) r ' • 2a 

2 j ^ 

Character* 

R 

C 

L 

L 

L 

c 

R 

R 

R 

Remark 


/3 = a 

k = 1 


f3 = 1 

k = 0 


a - 1 

k = 1 


/3 - a 

k = 0 



*The abbreviations are as follows: Revolution, Oitical, and Libration. 
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THE CASE OF MODERATE ECCENTRICITY 


In this case, we expand a and /3 into power series in j . From the expressions of Equation 9 it fol- 
fows, by neglecting high order terms, that 


- a Ail zo\ i 

T J \4 a 2 /3 0 

> 

n = • /lA1_ _Zl A 3 a 3 A 
J \4 a 2 /3 0 /3 0 a 2 2 a ^2 j 

Here, except for the critical case, we have only one of two cases: ib or ic 


Case ib: Libration 


(16) 


From Equations 8 and 10 we obtain 


du 


from which, expanding it into a power series in a in , we obtain 


3 “4 7 a 3 

8 2 J^Xo* 

2 


= 2a. fi 1/2 dr , 


(17) 


du 


3 a _4 15 a 3 

l 2 ” a 


1 ~M 1/2 Xo + + 


2a /a 1/2 dr , 


(18) 


up to the order ^ - 0( j ) , which corresponds to j 2 in Equation 1. We also have, from Equations 7 
and 8, 


Now, putting 


— a 1/2 v 

a„ M ao 


2^7^ 1/2 Xo + 1 


a 
2 a 


4 5 

+ 

2 


8 a 2/^X0 

a 2 


3^4 15 a 3 

ct 2 - " 2 a 2 + 8 a 2 ’ 

into Equation 18 and using Xo = cn\u to the required order, we have 

du ^1 + o-j fj. 1/2 cn Xu + <r 2 /n cn 2 = 2a 2 ,a 1/2 dr 

s du, . 


( 19 ) 
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It is known that (Reference 9, pp. 58 and 62) 


dv - k 1 sin (k sn v) , 


0 < k 2 < 1 


cnv 2 dv = k -2 [e(v) -k' 2 v 


where E(v) is the elliptic integral of the second kind such as 


E( 


- ■ f 


dn 2 v dv 


Accordingly, 


A.Uj - A (1 + cr 2 fisju + CTj k" 1 sin” 1 ( k sn Au ) + k~ 2 <r 2 fj.Z{ Au ) 


v , 


= 1 + k -2 (eK _ 1 - l) , 


E = 


/*tt/2 

(l - k 2 s in 2 v) 1/2 dv 
Jo 


77 / 11 , , 

2 F f 2 ’ 2 1 1 ! k 2 ) , 


and K is given by Equation 15. 
Equation 20 has a form 


Au„ = A (l + a 2 /as) u + (periodic terms in Au with the period 4K). 


Therefore, by inversion, it follows that 


Au - A 


Uj + ^periodic terms in A t 



where 


Thus 


X x - k (l +a 2 /xs ) -1 


Xu = vt - cr 1 k" 1 sin - 1 (k sn vtJ 

+ o-j 2 /x cn vt k -1 sin -1 (k sn vt ) - cr 2 /xk “ 2 Z( vt ) 


where 


( 22 ) 


V = 2 a J y£T Xj 

- 2 a // 2 j^ 2 r ^ 2 (a-/ 3 )l / 2 (l-<r 2 ( as) , 

andk 2 = ( 1 — > 8 ) (a: - >S) “ 1 . The other notations, a 1 and <j 2 , are given by Equations 19; s by Equation 21 
and a, /3 by Equations 16. Then, 


sn 2 (Xu, k) 
a “ P cn 2 (Xu,k) 

x I“3 fW 2 , 

cn Xu - 2 ^7 V ^7 cn Xu 

w (\ fix y o i a 4^o 5 a 3 2 ^o\ , ( 

+ ^2 8 ^ j cn ' 

Therefore, we have, neglecting higher orders, the following forms as the solutions of Equation 1, 



jf (l T To 
T a 2 fio /^o 


- 1^Y /2 

1 + Wo 

(1 T 2 

27o 

^2 3 

a 3^l\ 

cos “ " ^oj 

\2 a 2 /3 0 

O 

TO . 

[ 

" a 2 + 2 

a 2 2 y 


/XV /2 JL llL Zl h 3 a 3 Pl\ 3 

VV 2 !/3 0 \4 a 2 /3 0 /3„ a 2 + 2 ^2 ^ sn 


cn Xu - -9 — — i cn 2 Xu 
a 2 

_1 AH . 2 \ + /Jr jr !\ P 1 y o 

2 a 2 /3 0 J sn Xu y a 2 /?„ ^4 a 2 /3 Q £ 0 


1 

a 4 A> 

5 “3 2 /3o\ 

2 

*/ 

® a 2 2 a 2 y 


£ - 


*-£ 


jr /1 T 2 yp 
0 a 2 /^0 ^0 
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where Xu is given by Equation 22. The order of magnitude adopted here corresponds to the second 
order of Equation 7. 

Case ic: Revolution 

In this case a formula similar to the aforementioned one can be given: 

du (l + oq dn Xu + <j 2 /j. dn 2 Xu) = dUj , 

from which we have 

\uj = Xu + <7j am Xu + cr 2 /j.E(Xu) 

= ^1 +CTj m 172 f K' 1 +a 2 /xEK -1 j Xu (23) 

+ a-j /u 1/1 2 ^am Xu - 2K + a 2 fj. k~ 2 Z( Xu ) . 

Putting 

v - 2a 2 jj. 1/2 ( 1 - J3) 1/2 |l + o-j /j. 172 ~2 K~ 1 + a 2 /i EK” ^ “ 1 , 

we have, by inversion, 


Xu - vr - cr ; /x 17 2 I am vt — ^ K 


+ a 2 n |dn vt - -j K *^|am vr - y K' 1 vrj - cr 2 fx Z( vt ) , 


where 


k 2 = (a-/3)(l - /3) ~ 1 . 

We note that <j 1 /j . 172 am Xu, the second term of the right side of Equation 23, also contributes to the 
change of the period- contrary to Case ib, where the corresponding term, a 1 fi 1/2 k' 1 sin -1 (k snXu) , is 
purely periodic provided 0<k 2 <i. It is easily seen that, roughly speaking, Case c corresponds to 
Case b where k> l. At any rate, for the former case we have the following solutions of Equation 1: 


cos co - sinciq 


sn Xu 


1 + y j 

(i 


X 0 

^2 3 

a 3^i\ „ 

V 

a 2 /3 0 

"^0 

a 2 + 2 

, civ 1 Xu 

a 2 J J 
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p 


P + X 


d x 1 • 

dn Xu - 2 j — 


' Xu - 


r “3 

7 dn 2 Xu 

n ^ 


~ j 


Hr (i P? y 0 \ 2 

V * 2 y 4 a 2 / S 0 " / 3 jsn Xu dn Xu 




5 

+ 8 



Xu . 


GENERAL SCHEME OF THE FURTHER APPROXIMATIONS 

In order to solve the equations with higher order terms, it is first necessary to solve an equation 
such as Equation 6 which contains these terms. Even if the original equations of motion are limited 
to finite terms, the Equation 6-type of equation is, in general, an algebraic equation with cos 2 & as the 
unknown. When combined with the second equation of Equation 4, we find that (d/d r) (cos 2 o>)is equal to 
an algebraic function in cos 2 o>. The solution of a differential equation of this kind would involve 
Abelian integrals— many valued functions, more complicated than the elliptic integrals. Therefore, 
the solution becomes so complicated that we cannot see the general character of the motion at a first 
glance. However, if we restrict ourselves to the solution of a real function, the complicacy would be 
much reduced. 

In this section we are concerned with the cases where the character of motion is quite similar to 
the former cases, except for slight changes caused by the higher order terms. 

Now, we consider the equations of motion in the form: 

d X <9<P 

dT ~ da > 1 


where 


dt2J i 

dr 


dV 

d x 


4* = j a 0 (sin 2 aq , jj + Z! a n (sin 2 ^, x n 

and a n (n = 2,3, • •• } is a power series in j , a 2 ^sin 2 u> 1 , 0^ = a 2>0 and is a constant which corresponds 
to the previous case if a 2 0 - a 2 . 

From the integral V - jT = 0, we obtain 

4 a 2 2 x 2 (i + 2 a n* x n " 2 ) - 4 a 2 i( r_a o) = r 2 - 

from which 

2a 2 X = V + L b n v n • 
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Thus, we have 


(5) = dr 

or 

rf 1 ^1 + 2 c n 7j n j “ 1 da>j = constant • dr 

If the solution of v ~ l = du* could be obtained in a simple way, as in the previous cases, then 
using u* as an intermediary parameter would be a direct extension of the previous cases. But this 
integral is, in general, an Abelian integral as we have just stated. A general treatment is not suitable 
for practical purposes. As an alternative, we shall consider the following situations. 

a : Case of moderate eccentricity 

Let a Q (s in 2 a>j , 0) = a 0 , sin 2 o)j, a„ j > 0 J then the Solution Of a Q (sin 2 o)j , 0 ) -r = 0 is 

sin 2 ajj = (a 0tl ) _1 r . (24) 

Now, we denote one of the solutions by a -1 : 

a 0 (sinojj, j) - T = 0 , 

the first approximation of which is Equation 24. We then have 

4 a 2 j P - a 0 ^sin 2 oq, = fi 1 (l - a sin 2 aijj ^1 + 2 A n sin 2n a , 

where = 4Ta 2 0 j . Therefore, 

V = fFf £ (i +2 - fi n sin^tJi) , 

where E? - 1 - a sin 2 <n 1 and 

£ -1 1 + I/x” /2 f n (sin 2 £ l ’j“ 1 dc^ = /Xj 172 dr . 

Now 

E' 1 dco t = du . (25) 

Here we again divide the problem into two subcases: 

Case ib: a > 1 
We have 

= sn^Xu, a~ 1/2 j , 


yaTsin a> i 
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§ - y 1 - a sin 2 ojj 


then 


- cn (\u , a 1/2 ) , 
X = a 1 ^ 2 ; 


[l +S A 6 n/2 f n (a -1 sn 2 a 1/2 u) (cn n a 1/2 u)]' 1 du = /u, " 2 dr 


or 


[l + S M n/2 g n (a 


u) 


, 1/2 


-] 


du 


^ 2 dr , 


(26) 


where g n is a power series in j with the coefficients polynomials in a’ 1 sn 2 \u. Therefore, the in- 
tegration of Equation 26 is expressed as 


5 -f 


+ 2 B I cn 11 A.u du - 


1/2 , 


(27) 


where B n is a numerical constant. This integration is carried out by the use of the Jacobian elliptic 
function and the Jacobian zeta- function (see Reference 10, p. 164). Thus, the final form is expressed 
by 


Xu +. a periodic function of Xu = constant • r , 
where the constant factor in the right side is /x t 1/2 x[l + o( j )]. This is an extension of Equation 20. 
Case ic: a < 1 
From Equation 25 we have 

sin o)j = sn (u, a 1/2 ^ , 

£ - j/ 1 - a sin 2 aij = dn ^u, a 1/2 j ; 

therefore, 


[l +2 / a n/2 f n |sn 2 u) dn" uj' 1 du = dr 
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or 


|l + £ ;u n//2 g n (sn 2 u) dn n uj «du - / u. 1 1,/2 dr . 

In this case we have, instead of Equation 27, 

u + 2C n | dn n udu 3 , 

where the integration will be carried out, as above, by using the Jacobian elliptic function and the 
Jacobian zeta- function. Thus we have 

u + a periodic function of u = constant • t, 


where the constant factor in the right side again is 1/2 [l + 0( j )]. 

fi : Case of small eccentricity 

In this case the most important factor in r - a 0 (sin 2 is not the 1 - a sin 2 &> 1 term, as in the 

previous case, but is (l - a sin 2 oj^J (l - /3 sin 2 , where both a and P are of the order of 1. Therefore 
we cannot expand each of the above factors into power series beginning with a constant term. 

Or more concretely, hereafter a 01 is assumed to be of the order of j : a 01 = j a' 0 x and, also, it 
is assumed that 

lim j~ 1 a 0 (sin 2 a> r j\ = a/ x sin 2 a)j + a Q 2 sin 4 , 


where a 0 ' 4 corresponds to fi 0 ' , and a 0 2 to - (1/4) // a "/ Also, we denote by a -1 , p~ l the solutions of 


a 0 (sin 2 co v jj - j T' - 0 


which approach the solutions of 


T' - a/ j sin 2 c*J 1 - a Q 2 sin 4 aq = 0 
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As before, we have 


4a, j 


jr' -a 0 (sin 2 <Oj , j) = /ij - a sinoij) -/? sin 2 (l + 2 A n sin 2n , 


where ^ = 4r* a 2 0 j 2 . Therefore, 

= M! 1/2 X 0 (l +2^ n sin 2 " ojj^ , 

X 0 2 = (l-asin 2 d>^ (l -fi sin 2 

and, as before, 

X" 1 1+2/x" 72 f n (sin 2 co^j x 0 n J _1 Auj = /j.^ 2 dr . 

Since x 0 _1 = du, we have, for case i (except for the critical cases): 

sn 2 (Xu, k) 

x = sin 2 a) = ~ 77 TV 

1 p - q cn 2 ^Xu j k) 

where X, k , p, and q can be obtained in Table 1, according to the respective cases. Further, 

(p - q) cn 1 Xu dn^ 2 Xu 


(28) 


(29) 


p - q cn 2 (XUj k) 


(30) 


where e 1 , e 2 = 0 or 1, but both are not zero at the same time. 
Now, from Equation 28 we have 


[l + S/Xj" 72 f n (x)x 0 n j -1 du - Mi 1/2 dr 


or 


[l +2^ 1 n/2 g n (x) x 0 n ] du ' ^x' 2 dT 

where g is a polynomial of x . Thus, for n even, g n has the form 


(31) 


g n ■- G ! ( sn2 ;vu ) * (32) 

where G x denotes a rational function of the argument. The integration of Equation 32 produces an el- 
liptic integral of the third kind, in general, unlike Case a (Reference 10, p. 164). On the other hand, 
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for n odd, 


g n - G 2 [sn 2 A.u) cn 1 A.u dn 1 X.u , 

the integration of which requires only the Jacobian elliptic functions. 

For Case ii, however, the equations corresponding to Equations 29 and 30 have the different 
forms; but all the functions g n have a single form such as Equation 32. 

In any case, we have, after integration of Equations 28 or 32, the form 

\u + a periodic function of A.u = constant ■ r, 

from which we obtain 


= constant • r + a periodic function of constant ■ r. 

It should be noted that the odd power terms in Equation 31 do not contribute any change to the 
period of libration, where e t = 0, which corresponds to the case of libration. 

DISCUSSION 

The denominator of the right side of Equation 11 is a square root of a quartic of x. A general 
treatment of such a case was carried out by Andoyer (Reference 11) and extended by Hagihara (Refer- 
ence 12) in connection with the libration problem of asteroids (see also Reference 13). Both of them 
used the ^-function of Weierstrass; nevertheless, we feel the integration form described here is more 
suitable for this special case, owing to the fact that only real functions are involved. For example, if 
p - 0 in Equation 10, 


da)j 

yi - k 2 sin 2 ^ 1 


du , 


where k 2 = a, from which, if 0 < a < 1, we have sinu, = snu, which is, of course, a limiting case of 
ic, etc. Thus the present work represents some natural extensions of the derivation of well-known 
Jacobian elliptic functions. 

For the numerical computation, integrating Equation 17 in the literal form may not necessarily 
be the best way, for a direct numerical integration of Equation 17 is more effective, since its de- 
nominator generally does not become very small— otherwise, the numerical integration would lose its 
validity. However, the problem of inversion of a numerical function still would remain. 


25 



CONCLUSION 


Equations 1, which have terms up to second order of j, have been solved both for the case of 
small eccentricity and for moderate eccentricity. In the former, the intermediary solutions them- 
selves are sufficient for the present discussion, since the solutions have an accuracy of the order of 
1 for cos « (neglecting the order of j) and of the order of j for P (neglecting j*). This accuracy cor- 
responds to the order of j 2 in the Hamiltonian, since p has the factor of j and the argument of the 
solution is neither r alone nor multiplied by a constant factor of order of 1, but r multiplied by that 
of order j . 

On the other hand, in the second or normal cases, the solutions have the accuracy of order j for 
cos co and j 3/ 2 for p, which corresponds to j 2 in the Hamiltonian. 
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Appendix A 


C d’ C n’ d n 

e f , e g , e h 
E 

f i' 8j- h i 

Jn 

k 

K 


t 

T 


a, / 3 , a. , / 3 . , 7 . 

r = jr- 


p 

O) 


aij = it/ 2 - cd 




List of Symbols 

Jacobian elliptic functions (see pages 12, 8, and 9) 

Parameter of Weierstrassian elliptic function (see page 6) 

Complete elliptic integral of the second order (see page 17) 

First order elliptic function associated with the function (see Reference 6 
and page 7) 

Harmonic coefficient of the earth's potential (see page 2) 

Modulus of the Jacobian elliptic function (see pages 7 and 11) 

Complete elliptic integral of the first order (see page 14) 

Weierstrassian elliptic function (see page 6) 

Jacobian elliptic function (see page 8) 

Time 

Period of the solutions (see page 14) 

u will equal t, except for the constant factor, if higher orders are neglected 
(see page 13) 

Constants 

Energy constant (see page 4) 

Constant factor in the argument of the Jacobian elliptic functions (see 
page 12) 

Variable related to the eccentricity (see Reference 5 and page 2) 

Mean argument of the perigee 

Semi-periods associated with the Weierstrassian elliptic function (see page 7) 
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